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LONG-RANGE,  DEEP-OCEAN  PROPAGATION  OF 
15-Hz  CW  ACOUSTIC  SIGNALS 


DESCRIPTION  OF  THE  EXPERIMENT 


A ship  towed  two  continuous-wave  (CW)  sources  operating  near  15  Hz  at  140-m  depth 
and  114  Hz  at  21-m  depth,  dropped  explosive  charges  on  a regular  schedule,  and  made 
environmental  measurements.  The  acoustic  signals  were  received  on  an  omnidirectional 
sound-channel  hydrophone  which  is  suspended  up  firom  the  bottom  in  deep  water  near 
Midway  Island.  The  exercise  lasted  14  days  and  was  divided  into  five  phases.  During  the 
exercise  the  ship  made  47  expendable  bathythermograph  (XBT)  casts  and  kept  a contin- 
uous bathymetric  record  by  using  a precision  fathometer  (PFR). 


Figure  1 is  an  abstract  of  the  experimental  geometry.  The  first  phase  was  a great  circle 
transit  from  the  hydrophone  area  to  a range  of  700  km  made  to  determine  the  environmen- 
tal parameters  along  this  track.  Throughout  the  remaining  phases,  the  acoustic  sources  were 
operated  and  the  ship  attempted  to  maintain  a constant  7 knot  speed.  In  the  second  phase, 
the  ship  continued  along  the  extension  of  the  great  circle  track  of  phase  one  to  a range  of 
1700  Im.  In  the  third  phase  the  ship  reversed  course  and  returned  to  a range  of  700  km. 
The  fourth  phase  repeated  the  second  phase.  The  intention  in  the  fifth  phase  was  to  have 
the  ship  steam  in  the  direction  transverse  to  the  line  of  acoustic  propagation.  In  fact,  the 
ship  proceeded  roughly  along  a 600-km-long  straight  line  track  from  the  terminus  of  phase 
four  at  a 1650-km  range  to  a point  1550  km  from  the  hydrophone. 


Fig.  1 — Abstract  of  the  ship’s  track 
labeled  to  show  the  five  phases  of  the 
experiment  and  the  nominal  range 
from  the  receiver  at  the  beginning  of 
each  phase 


receiver 


The  CW  sources  were  controlled  by  a crystal  standard.  At  the  receiver,  each  CW  signal 
was  band-pass  Rltered,  heterodyned  to  low  frequency,  low-pass  filtered,  and  sampled  once 
every  two  seconds.  Signals  from  the  explosives  were  analog  recorded. 


A description  and  a listing  of  the  environmental  data  acquired  have  been  published 
[1] , as  have  the  results  from  the  shallow-explosive  signals  [2,3] . The  15-Hz  CW  signals  from 
phase  five  have  been  analyzed  and  interpret^  in  terms  of  signid  coherence  along  a line 
transverse  to  the  direction  of  propagation  [4] . This  paper  reports  the  15-Hz  CW  trans- 
mission loss  data  obtained  in  phase  two,  three,  and  four. 


Manuscript  rabmitted  December  29, 1978. 
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ENVIRONMENTAL  DATA 
Navigation 

The  basic  navigation  data  consisted  of  Navy  Navigation  Satellite  System  fixes  as 
available  plus  Loran  C fixes  every  15  min.  Some  information  about  ship’s  speed  was  used  to 
improve  the  navigation.  We  are  confident  that  the  ship’s  range  is  accurate  continuously  to 
±1  km  during  phases  two,  three,  and  four  which  are  the  subject  of  this  paper.  This  accuracy 
is  confirmed  by  the  range-sensitive  analysis  [2]  performed  on  the  explosive  signals.  Figure  2 
shows  details  of  the  navigation  for  phases  two,  three,  and  four. 


Bathymetry 

The  bathymetry  determined  from  the  PFR  records  for  an  assumed  sound  speed  of 
1500  m/s  is  plotted  in  Fig.  3.  Although  not  shown  in  this  figure,  the  bottom  rises  in  the 
vicinity  of  the  hydrophone  to  a depth  of  about  3000  m.  A brief  survey  of  the  seamount  at 
400  km  indicated  that  it  rose  to  a minimum  depth  of  3912  m just  south  of  our  track.  The 
feature  at  1530  km  showing  a marked  depression  on  both  sides  of  the  raised  central  section 
belongs  to  the  Mendocino  Fracture  Zone. 


Thermal  Structure  of  the  Water  and  Sound  Speed  Profiles 

Two  types  of  Sippican  XBT  were  employed  — a mod  T7  having  a design  depth  of  760  m 
and  a mod  T5  having  a design  depth  of  1800  m.  The  intention  was  to  sample  the  more  stable 
deep  layers  less  frequently  than  the  unstable  shallow  layers  by  interspersing  the  more 
expensive  mod  T5  among  the  mod  T7  XBT  casts.  The  results  are  disappointing  in  that, 
although  either  the  T5  or  T7  casts  taken  as  a group  are  consistent,  the  two  groups  are  not 
consistent  with  one  another.  The  T7  casts  resemble  the  T5  casts  down  to  about  400  m,  but 
the  T7  casts  indicate  a significantly  lower  temperature  from  400  m to  760  m.  Based  on  the 
facts  that  the  T5  casts  agree  better  with  archival  data  [1]  (the  very  low  temperatures 
indicated  by  the  T7  casts  from  400  m to  760  m disagree  with  archival  data)  and  that  the 
pattern  of  explosive-shot  arrivals  was  successfully  modeled  by  using  T5  profiles,  the  T7 
casts  have  been  set  aside.  The  bathymetric  measurements  and  the  sound  speed  profiles 
derived  from  T5  XBT  casts  along  the  1700-km  radial  track  of  this  experiment  are  summa- 
rized in  Fig.  4. 


ACOUSTIC  DATA 

Figure  5 contains  a representative  sample  of  the  heterodyned  and  digitized  signal 
received  from  the  15-Hz  source.  The  Hilbert  transform  of  this  digital  signal  was  computed 
and  used  in  conjunction  with  the  navigation  to  obtain  the  signal  amplitude  at  0.5-km 
intervals.  This  spatial  sampling  rate  exceeds  the  Nyquist  rate.  The  amplitude  is  plotted  in 
Fig.  6 as  transmission  loss  vs  range.Figures  7 and  8 are  smoothed  versions  of  Fig.  6 obtained 
by  plotting  the  average  received  intensity  over  7-km  and  35-km  intervals,  respectively,  as 


- Measured  bathymetry  for  an  assumed  sound  speed  of  1500  m/s  in  (a)  phase  one,  (b)  phase  two,  (c)  phase  three,  and  (d)  phase  four 
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Fig.  4 — Environmental  summary  for  phases  one,  two.  three,  and  four. 
The  SOFAR  axis  and  critical  depth  are  shown  as  dashed  lines. 


transmission  loss  vs  range  at  the  center  of  the  interval.  For  the  frequency,  ranges,  and  nine- 
day  time  interval  of  the  experiment.  Fig.  8 gives  an  indication  of  the  temporal  stability 
of  the  mean  transmission  loss  curve  defined  by  a 35-km  range  average  of  intensity. 


Fig.  5 — A typical  60-min  sample  of  the  digital  acoustic  data  from  the  15-Hz  source  in  phase  three,  llie 
digital  signal  was  obtained  by  filtering  and  heterodyning  the  received  signal  and  then  sampling  once  every 
2 s.  The  high  signal-to-noise  ratio  is  evident  as  is  the  stability  of  the  phase  and  amplitude  as  the 
ship  proceeds  at  7 knots  in  the  direction  of  the  hydrophone. 


4 


ISS  IdB) 


NRL  REPORT  8302 


V ':  : I .'  ‘fil  l\i  i*  .'*■ 

i i If  T,!  '!  |.!  ' lii' 

■i  ! ‘ . i!  I Mil  ri 


RANGE  Iktn  x 100) 

Fig.  6 — Transmission  loss  vs  range  of  the  15-Hz  signal  for  (a)  phase  two,  (b)  phase  three 
and  (c)  phase  four.  The  spatial  sampling  rate  is  two  samples/km. 
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A noteworthy  feature  of  Fig.  7 is  the  presence  of  zones  of  constructive  interference 
at  53.5-km  intervals  with  weaker  such  peaks  appearing  halfway  between  them.  Figure  7(a) 
shows  these  features  particularly  well  in  the  interval  between  1000  km  and  1600  km.  The 
fact  that  the  ship  took  48  h to  traverse  this  interval  indicates  a measure  of  long-term  stabil- 
ity not  only  in  the  broad-area  aveif^^e  transmission  loss  but  also  in  some  of  the  details  of 
the  interference  pattern.  The  time  interval  from  the  beginning  of  phase  two  to  the  end  of 
phase  four  includes  nine  full  days;  yet,  some  features  of  the  interference  pattern  are 
repeated  in  all  three  phases. 

In  order  to  understand  the  observed  periodicity,  we  employed  a normal  mode  model 
of  the  acoustic  field  which  required  that  the  medium  be  specified  by  a single  sound  speed 
profile.  Because  Fig.  4 indicates  a significant  change  in  the  ocean  sound  speed  structure 
at  a range  of  about  1000  km,  two  normal  mode  calculations  were  performed.  The  first 
calculation  used  the  single  sound  speed  profile  of  Fig.  9(a)  which  summarizes  the  water 
properties  from  0-  to  1000-km  range.  The  second  calculation  used  the  profile  of  Fig.  9(b) 
which  summarizes  the  water  properties  at  ranges  in  excess  of  1000  km.  In  both  cases,  the 
bottom  was  taken  to  be  a flat,  rigid  reflector  at  a depth  of  6000  m.  The  surface  was  mod- 
eled as  a perfectly  reflecting  pressure  release  boundary,  and  no  volume  attenuation  was 
considered.  The  neglected  effect  of  volume  absorption  amounts  to  about  0.5  dB  at  1000  km 
for  15  Hz  (5) . To  describe  the  modes  and  their  properties  we  employ  the  ray -mode  analogy 
and  number  the  modes  in  the  conventional  manner  so  that  mode  1 corresponds  to  the 
SOFAR  axial  ray.  For  both  profiles.  Fig.  9,  (a)  and  (b),  it  is  found  that  modes  1-19 
correspond  to  RR  rays,  modes  20-27  correspond  to  RSR  rays,  and  modes  numbered  greater 
than  27  correspond  to  BRSR  rays. 

A brief  digression  is  in  order  at  this  point  to  consider  the  region  of  validity  of  this 
model.  If  the  observed  acoustic  field  depended  strongly  on  the  detailed  acoustic  properties 
and  structure  of  the  ocean  bottom,  the  model’s  validity  would  be  in  grave  doubt.  As  we 
shall  see,  the  observed  acoustic  field  is  composed  primarily  of  trapped  modes  which  are 
relatively  independent  of  the  bottom  properties.  The  higher  order  bottom -sensitive  modes 
excited  by  the  source  are  removed  from  the  propagating  acoustic  field  by  their  interaction 
with  the  ocean  bottom.  The  normal  mode  calculation  is  valid  in  the  sense  that  the  calculated 
trapped  modes  model  the  experimentally  observed  field.  Higher  order  bottom -sensitive 
modes  are  calculated  in  the  model  to  demonstrate  that  if  something  more  than  just  the 
trapped  modes  were  present  in  the  propagating  field,  the  properties  of  the  calculated  field 
would  disagree  with  those  of  the  measured  field.  Simplicity  and  ease  of  calculation  were 
the  reasons  for  selecting  the  particular  bottom  model  used  to  calculate  higher  order  modes. 

In  the  following  paragraphs  various  comparisons  are  made  between  the  experimental 
data  and  the  calculated  acoustic  field.  We  consider  first,  the  average  received  intensity  over 
a 35-km  range  interval;  second,  the  average  received  intensity  over  a 7-km  range  interval; 
and  third,  the  spectrum  of  the  range-dependent  intensity  corrected  for  cylindrical  spreading. 

Figure  10  shows  a comparison  between  the  35-km  average  intensity  of  phase  two 
(Fig.  8(a))  and  the  calculated  average  for  the  profile  of  Fig.  9(a).  It  is  clear  that  the  RR 
modes  alone  are  capable  of  accounting  for  the  average  received  intensity.  If  we  consider  the 
possible  errors  in  source  level  and  system  response,  however,  this  comparison  between  model 
and  data  is  taken  to  be  more  suggestive  than  conclusive.  Figure  11  shows  a similar 
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Fig.  9 — Sound  speed  profiles  used  in  the  normal 
mode  calculations.  Profiles  (a)  and  (b)  charac- 
terize the  environmental  properties  for  ranges 
less  than  and  greater  than  1000  km,  respectively. 


comparison  between  the  data  and  the  calculated  average  intensity  for  the  profile  of  Fig.  9(b) 
Again  the  RR  modes  are  sufficient  to  explain  the  received  intensity. 

Figure  12  shows  the  average  intensity  received  over  a 7-km  interval  as  calculated  for 
the  profiles  of  Fig.  9.  In  the  range  interval  1300  to  1700  km.  Fig.  12  shows  zones  of 
constructive  interference  at  53.5-km  intervals  with  secondary  peaks  midway  between.  The 
RR  modes  producing  these  peaks  have  deep  turnaround  points  between  3000-  and  3600-m 
depths.  Such  modes  would  propagate  with  little  interference  from  the  topography  of  Figs.  3 
and  4.  Although  Fig.  12  is  conditioned  by  such  caveats  as  (i)  only  modes  1 to  12  were 
summed  to  obtain  the  pressure  (which  causes  the  calculated  transmission  loss  to  be  higher 
than  that  of  the  experimental  data)  and  (ii)  the  ranges  of  the  calculated  and  experimental 
transmission  loss  peaks  do  not  match,  it  is  significant  that  the  observed  periodicity  in  the 
experimental  transmission  loss  can  be  modeled  by  the  same  subset  of  RR  modes  for  either 
profile  in  Fig.  9. 
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The  normal  mode  model  of  acoustic  intensity,  <p^>,  when  only  a finite  number  N 
of  modes  propagates  between  a source  and  receiver  at  fixed  depths,  is  given  by 

<p2>  = r-l  ^ A„A^co8(k^-k^)r,  (1) 

fn,n*  1 

where  describes  the  excitation  of  mode  m and  r is  the  range.  The  radial  wave  number 
is  a monotonically  decreasing  function  of  the  mode  number  m.  Spectrum  analysis  of 
yields  information  regarding  the  distribution  of  components  present  in  the 

data.  The  largest  difference,  corresponds  to  interference  between  the  excited, 

propagating  modes  having  the  lowest  and  highest  mode  numbers.  Thus  another  estimate  of 
the  number  of  propagating  modes  present  in  the  data  can  be  obtained  by  a comparison  of 
the  falloff  in  the  data  spectra  of  <p^>  ask„  - k„  increases  with  the  falloffs  in  the  spectra  of 
<p^>  calculated  for  the  profiles  of  Fig.  9 for  various  numbers  of  propagating  modes. 

Figure  13  shows  spectra  of  r<p^>  for  intervals  of  512  km  beginning  at  the  minimum 
range  in  each  phase.  Figures  14-17  show  calculated  spectra  of  r<p2>  for  the  profile  of 
Fig.  9(a).  These  spectra  were  obtained  from  the  calculated  values  of  <p^>  at  0.5-km 
intervals  from  700-  to  (700  + 512)-km  ranges.  All  spectrum  levels  are  given  in  decibels 
relative  to  the  same  arbitrary  level.  If  we  compare  Fig.  13  with  Figs.  14-17,  it  appears  that 
the  highest  mode  number  of  a propagating  mode  is  about  19.  Again  the  conclusion  is  that 
no  bottom-interacting  modes  appear  in  the  experimental  spectrum. 

This  type  of  analysis  is  insensitive  to  errors  in  source  level  and  system  response,  pro- 
vided that  the  errors  do  not  change  with  time.  During  the  course  of  the  experiment  and  in 
the  processing  of  the  data,  no  changes  were  made  in  the  responses  of  the  low-frequency 
receiving  and  processing  systems.  The  15-Hz  source  level  could  not  vary,  because  the  source 
has  fixed-displacement  pistons.  The  spectrum  analysis  is  conservative  in  determining  an 
upper  bound  on  propagating  mode  number  in  that  it  assumes  temporal  stability  in  the 
transmission  loss  curve  over  the  approximately  40  h required  for  the  ship  to  go  512  km. 
Temporal  variations  in  transmission  loss  with  periods  of  10  min  or  less  (corresponding  to 
range  intervals  of  2 km  or  less  for  a ship  at  7 knots)  would  appear  in  the  data  spectrum  in- 
distinguishable from  contributions  of  BRSR  modes.  In  addition,  the  spectrum  analysis  is 
conservative  in  estimating  an  upper  bound  on  propagating  mode  number,  because  we  have 
compared  the  data  spectra  with  spectra  calculated  by  using  a normal  mode  model  in  a 
range-independent  medium.  Any  range  dependence  present  in  the  ocean,  whether  in  the 
sound  speed  profile  or  in  the  bathymetry  (cf.  Fig.  4),  produces  coupling  among  the  modes 
which  form  a mathematical  basis  for  the  acoustic  field.  Such  an  increase  in  the  number  of 
basis  functions  would  appear  in  the  spectra  indistinguishable  from  contributions  of  higher 
order  propagating  BRSR  paths. 

Similar  conclusions  can  be  drawn  firom  Fig.  13  by  comparing  it  with  an  even  simpler 
model  of  the  acoustic  field.  Suppose  there  are  N propagating  modes  with 

“ 1 (2) 
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Fig.  14  — Power  ipectral  density  of  r<p*>  for  ranges  700  to  (700  512)  km  and  for  the  first  12  modes  of  the  profile  of 

Fig.  9(a).  TTie  appearance  of  spectral  components  to  the  right  of  k„  “ “ 2n/2307  m“*  in  this  figure  is  due  to  leakage 

associated  with  the  use  of  a rectangular  data  window  and  to  finite  computer  word  size. 
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and 


for 


♦ ma 

in 


(3) 


"»  = 1.2 N.  (4) 

where  k and  a are  constants.  Equation  (1)  becomes 
N 

r<p^>  • 2 co${m-n)ar 

(5) 


• N * 2{N- 1)  cosar*  2(N- 2)  cos  2 a r * ...  * 2 cos  (N-l)  a r, 

which  is  in  the  form  of  a Fourier  series.  If  we  neglect  the  dc  peak,  the  average  power  spec- 
trum of  r<p2  > will  fall  off  by  15  dB  when 


or  when 


15  - 20  log 


2(^-1) 

2(N-j) 


~-l  - 10-»5/2O  (l-AT-l) 
/V 

If  we  assume,  conservatively,  that 


then  Eq.  (7)  implies 


Af  5, 


(6) 

(7) 

(8) 


j/N  > 0.86  (9) 

at  the  point  in  tlie  spectrum  where  the  average  is  down  by  15  dB. 

From  Fig.  13  it  appears  that  the  power  spectrum  of  r<p2>  is  down  15  dB  certainly  by 


2ir 

4000 


m 


-1 


(10) 


If  2a/4000  m'^  is  86%  of  the  value  of  which  corresponds  to  the  interference 

between  the  lowest  and  highest  numbered  propagating  modes,  then 


* 2ir/(4000  X 0.86). 


(11) 
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The  lowest  numbered,  excited,  propagating  mode  is  reasonably  approximated  (and  the 
approximation  has  been  validated  by  the  normal  mode  calculations)  as  the  mode  corre- 
sponding to  the  ray  which  turns  over  at  the  source  depth,  142  m.  The  average  speed  of 
sound  at  the  source  depth  was  about  1502  m/s.  This  gives  an  estimate  of  the  latest  wave 
number  as 


*m<w  = (12) 

Solving  Eq.  (11)  for  k^j„  and  using  Eq.  (12)  gives 

*min  * w/1502  - 2irf(4000  X 0.86).  (13) 

The  ray  corresponding  to  this  mode  turns  around  at  depths  where  the  sound  speed  is  given 
by 


=“1548  m/s 

for  our  source  frequency.  For  the  sound  speed  profiles  of  Fig.  4,  this  corresponds  to  an  RSR 
ray  which  turns  under  at  a depth  of  about  5400  m.  Figure  4 shows  that  the  bottom  depth 
exceeds  5400  m over  most  of  the  track.  This  simple  analysis  again  indicates  the  absence  of 
propagating  bottom-interacting  paths. 
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